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I will discuss the current status of the determination of ct(mz) and 
ct s (m z ). 

m 

^ INTRODUCTION 

a 

I ^ have been asked by the organizers of this conference to present a talk on 
the status of the electroweak sector of the Standard Model (SM). But frankly, 
4f"by the word 'status' one means 'how well theory and experiment agree with 
Ssach other', there is not much that I can tell you at the moment but that 
the SM works extremely well. Though the experiments at LEP/SLD and 
ewhere have pushed the experimental accuracy of many observables to a 
ere fraction of a percent, no significant deviation from the SM is yet to be 
ected. (I will discuss whether the 2.4(7 deviation in Ri, seen at LEP is 
nificant' or not later in this talk.) A detailed comparison between the SM 
igrr^d current precision electroweak measurements can be found in the talk by 

Dorothee Schaile in this proceedings (1). 
""^1 will therefore focus on the status of the determination of the values of 
a(mz) and a s (mz) instead. These quantities are used as inputs to calculate 
the SM predictions for the electroweak observables at LEP/SLD and any 
Certainty in them limits the accuracy of the theoretical predictions even 
ien the top and Higgs masses are fixed to given values. Consequently, any 
conclusion we may reach about the 'status' of the electroweak sector of the 
SM is actually contingent upon how well known these quantities are. 
Currently, LEP uses the following values for a(mz) and a s (mz): 

a(m z ) = 1/(128. 87 ± 0.12), 
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a s {m z ) = 0.123 ± 0.006. (1) 

The value of a(mz) is that determined from the <r(e + e~ — > hadrons) data 
from various experiments in Ref. (7), and the value of a s (mz) is that deter- 
mined from the hadronic event shapes, jet rates, and energy-energy correla- 
tion at LEP in Ref. (2). They are therefore independent of the Z line shape 
and asymmetry measurements at LEP and SLD. 

Note that the uncertainty in a(mz) arises from the incalculability of the 
hadronic contribution to the photon vacuum polarization which forces us to 
rely on the <r(e + e~ — > hadrons) data, and the uncertainty in a s (mz) arises 
mainly from our limited understanding of parton hadronization and our igno- 
rance of higher order QCD corrections. They are both manifestations of our 
limited ability to handle QCD. 

As mentioned above, these uncertainties in a(mz) and a s (mz) propagate 
into the SM predictions that are derived from them. For instance, for the 
SM with m t = 180GeV, m# = 300GeV, and Eq. 1 as input, the program 
ZFITTER 4.9 (3) gives the following predictions: 



sin 2 O 1 ^ = 0.2319 ± 0.0003 ± 0.0000, 

Rt = 20.766 ± 0.006 ± 0.040, (2) 

where R^ = T(Z — > hadrons) /Y{Z — > l + l~). The first error is that due to the 
uncertainty in a(mz) and the second error is that due to the uncertainty in 
a s (mz). The error in sin 2 due to the uncertainty in a s (mz) is negligible 
because QCD corrections only enter at 0(aa s ). 

On the other hand, the latest measurements of these quantities at LEP/SLD 
are given by (1) 

sin 2 O 1 ^ = 0.2315 ± 0.0004 

R t = 20.800 ± 0.035. (3) 

Comparison with the theoretical predictions Eq. 2 shows not only that the 
SM works very well, but also just how 'precise' the precision electroweak 
measurements at LEP/SLD have become. The experimental error is already 
comparable to or even smaller than the theoretical error 2 . With the experi- 
mental precision improving incrementally every year, clearly an improvement 
on the theoretical error is called for if one wishes to test the SM at an even 
higher level of accuracy. 

In the following sections, I will review the current status of the determina- 
tion of a(mz) and a s (mz) and discuss how the values quoted in Eq. 1 can be 
expected to improve or change. 



i e 1 Burkhardt et al. (1989) 
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' e ' Martin & Zeppenfeld (1994) 
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' s ' Swartz (1995) 
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a _1 (m z ) 

FIG. 1. Evaluations of a(mz) by different authors. All values have been rescaled 
to mz = 91.1887 GeV and the top quark contribution have been removed. 



THE VALUE OF a(m z ) 
Current status 

Quite recently, several authors have independently made attempts to reeval- 
uate the value of a(mz) through a careful reanalysis of existing <r(e + e~ — > 
hadrons) data (8-11). In comparison to the previous determination in 
Ref. (7), the hope was that the uncertainty in a(mz) could be reduced since 
some new data on <r(e + e~ — > hadrons) have been released (12-14), and the 
O(Oj) QCD correction to the cross section (15) together with a better de- 
termination of Q s (mz) (16) was now available. The resulting new values of 
a(mz) are shown in Fig. 1 together with a couple of older values. 

As you can see from Fig. 1, instead of finding a significant decrease in the 
uncertainty of a(mz), different authors found central values of a(mz) which 
differed from each other by as much as 2<7. Please note that all these values 
were obtained by applying different analyses on essentially the same set of 
data so the statistical errors are almost completely correlated. Therefore, 
large differences in the central values must be manifestations of systematic 
errors which have been underestimated. 

At the time of this conference, the two latest evaluations by Eidelman and 
Jegelehner (10) and Swartz (11) were not yet available so the disagreement 
among authors was rather pronounced. With the two latest numbers, the 
value of a(mz) shows signs of converging to a common value. 

In the following subsections, I will discuss how the analyses of different 
authors differ from each other so that we may get an idea of where the dis- 
agreements shown in Fig. 1 are coming from. 



2 Ref. (4) argues that the experimental error on Ri may be greatly underestimated. 
I will assume that that is not the case in the following. 
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The Definition of a(mz) 

The quantity whose value is usually quoted as that of the "effective QED 
coupling constant at the Z mass scale" is defined as 

a- 1 (mz) = a" 1 [l-A«M (4) 

with 

Aa(s) = 47raRe [W QQ (s) - W QQ (0)] , (5) 

where IL'qq(s) is the photon vacuum polarization function with only the light 
fermion contributions included. It has been customary to exclude the top 
and W contributions from Aa(mz) since the top mass was unknown until 
quite recently (17) (though care is needed when comparing results since re- 
cent authors include it) and the W contribution is also excluded to keep the 
definition of a(mz) gauge independent. These contributions are also numeri- 
cally small compared to the light fermion contribution since the W + W~ and 
ti thresholds are above the Z mass so that they do not contribute logarithms 
to the running of a(s) between s = and s = m 2 z . 



Contribution of the light fermions 



The contribution of the leptons to Aa can be calculated accurately in per- 
turbation theory and one finds 



&an eptons {m 2 z ) = V" — 

= y - 

^ 3-77 



--+$--&(3-$)ln 



Pi 



1+Pi 



3tt 
0.03142, 



In 



5 „ ( m 
- + O ^ 
3 \m', 



(6) 



where = ^/l — 4mj/m z . 

On the other hand, the contribution of the five light quarks (u, d, s, c, b) to 
Aa cannot be calculated perturbatively. Instead, unitarity and the analyticity 
of IL'qq(s) is used to write 



Aa 



(5) 

hadrons 



37T J 4 



S'(s-S')' 



where 3 



R(s) 



<r(e^ 



7* — > hadrons) 
-> 7* -> /X+/X") 



-127rlmn' {s) 



(7) 



(8) 



3 Note that the cross section in the denominator of the definition of R(s) is the 
cross section of electron pairs annihilating into massless muon pairs. 
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and the functional form of R(s) is extracted from experiment. It is this 
reliance on the experimental values of R(s), which are always accompanied 
by experimental errors, that we end up with a relatively large error on a(mz) 
even though the find structure constant a is known to extreme accuracy. 

The major problem associated with using the experimental values of R(s) is 
that the data points are only available for discrete, scattered values of s while 
one needs R(s) for all values of s to calculate AafJ ir(ms (s). Two methods 
have been used in the literature to deal with this problem. The first is to 
connect the data points directly with straight lines and perform trapezoidal 
integration (5,7,10), and the second is to guess the functional form of R(s) 
and fit it to the data (5,8,9,11). 

Both methods have their pros and cons. Trapezoidal integration is free of 
human prejudice about the functional form of R(s) but it does not take into 
account the experimental errors properly: sparsely distributed precise data 
points may not get the appropriate weight relative to the densely spaced data 
points with larger errors. Connecting two data points that are far apart with 
a straight line will also introduce errors. 

On the other hand, fitting a guessed functional form to the data has the 
advantage that experimental errors are easier to take into account. (Though 
care is needed in treating normalization errors (18).) However, the result 
will depend on the choice of the fit function and its parameterization, and 
systematic errors and biases will be introduced that are difficult to estimate. 

The two methods are often combined (e.g. fitting Breit-Wigner forms to 
the narrow resonances and using trapezoidal integration for the continuum) 
and are supplemented by the use of perturbative QCD for the high energy 
tail of R(s). 

The analysis of Burkhardt et al. 

At the time when LEP started running in 1989, the most accurate deter- 
mination of ^tt[aiire»i( m z) was that given by Burkhardt et al. in Ref. (5): 

*°kL m .{™l) = 0-0286 ± 0.0009. (9) 

(Actually, Ref. (5) reports the value of Aa[ 5 j rMJ (s) at ~Js = 92GeV. Rescal- 
ing to s/s = mz = 91.1887GeV gives the above value (6).) This value was 
calculated using the following three methods of integration and all three were 
found to agree with each other: 

1. Trapezoidal integration for the continuum and the p. 
Breit-Wigner forms for the narrow resonances. 

2. Trapezoidal integration for the continuum after a partial smoothing out 
of the R(s) data. (Details are not given in Ref. (5).) 
Breit-Wigner forms for the narrow resonances. 



6 




4 6 8 10 

12-94 W (GeV) 7854A1 

FIG. 2. Comparison of MARK I (black circles) and Crystal Ball (while diamonds) 
data. The solid line shows the perturbative QCD result. 



3. Breit-Wigner forms for the p and narrow resonances. 

Linear interpolation in */s for every few points in the continuum. 

In all three cases, perturbative QCD (at 0(a 2 s )) was used for R(s) above 
40GeV. Eq. 9 corresponds to 

a- 1 (mz) = a" 1 [l- Aa, eptons {m 2 z ) - Aa^^ml)] 

= 128.81 ±0.12 " (10) 

and results in an error of about 0.1% in a(mz). The region below the bottom 
threshold contributed about 1/3 of AafJ ir(ms (s)i and 80% of the error. 



The analysis of Jegerlehner 



The result of Ref . (5) was subsequently updated in 1991 by Jegerlehner (7), 
who was one of the original authors, in which the data from the MARK I 
collaboration (21) in the energy region 5-7 GeV were replaced by the more 
accurate data from the Crystal Ball collaboration (22). Fig. 2 shows the 
data of both collaborations in the energy range in question, together with the 
perturbative QCD result. The method of integration was the same as method 
No. 1 of Ref. (5) and the result was 

^hLonsiml) = 0.0282 ± 0.0009, (11) 
with no change in the size of the error. This corresponds to 

a -1 (m z ) = 128.87± 0.12 (12) 
which has been the standard value for the past few years. 
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FIG. 3. Fit to the continuum part of R(s) with (solid line) and without (broken 
line) correlations from normalization errors. The difference between the two in the 
3.6-5 GeV region is due to a different treatment of the higher resonances. 



The analysis of Swartz 

The recent analysis by Swartz in Ref. (8) fits a smooth function described 
by polynomials in W = s/s to the continuum part of R(s). Resonances were 
described with breit-Wigner forms as usual, and perturbative QCD (at 0(a^)) 
was used above 15 GeV. 

The major difference from method No. 3 of Ref. (5), or any other analysis 
that had used fitting functions, was that the correlations between data points 
within the same experiment due to the overall normalization error had been 
taken into account. The effect of this was that the fit produced a curve which 
was somewhat lower than when the correlations were neglected. See Fig. 3 
for a comparison. 

With a lower curve for R(s), the value Aa^ drons (m 2 z ) turned out to be 
smaller, and a~ 1 (mz) larger than previous analyses. Subtracting out the top 
quark contribution from the value quoted in Ref. (8), we obtain 

Aa£ 5) , (ml) = 0.02672 ± 0.00075, 

hadrons \ 6 1 ' 

and 

a _1 (m z ) = 129.07 ± 0.10, (14) 

which differs from Eq. 12 by 0.20 « 2cr. 

However, this analysis has been criticized (10) on the grounds that including 
normalization errors in the correlation matrix will introduce a bias when the 
normalization error and the number of data points are large (18). Swartz 
has subsequently updated his analysis to correct for this problem and also to 



(13) 
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FIG. 4. The new fit by Swartz. 



50 100 



include some new data which was missing from his first analysis (11). The 
new fit to the continuum part of R(s) is shown in Fig. 4, and results in 



Aa 



and 



( 5) , (ml) = 0.02752 ± 0.00046, 

hadrons \ Z, J ' 

aT x (m z ) = 128.96 ± 0.06. 



(15) 
(16) 



The difference with Eq. 12 has been reduced by 1/2, and the uncertainty has 
also decreased due to the inclusion of new precise data. 



The analysis of Martin and Zeppenfeld 

The analysis of Martin and Zeppenfeld in Ref. (9) distinguishes itself from 
all the other analyses in its extensive use of perturbative QCD. In table 1 I list 
the energy ranges in which different authors have applied perturbative QCD. 
(Though I also list the values of a s (mz) that were used, at the current level 
of accuracy, the difference is insignificant. Changing the value of a s (mz) has 
little effect on the resulting value of a(mz).) 

In the two energy regions 3-3.9 GeV and 6.5-oo GeV, Martin and Zep- 
penfeld express R(s) as the perturbative QCD value plus the J/^, ^, and T 
resonances. The DASP (19), PLUTO (20), MARK I (21), and Crystal Ball 
(22) data are all rescaled to fit the perturbative QCD result in these energy 
ranges, and the rescaled data is used to resolve a couple of resonances in the 
energy region between 3.9 and 6.5 GeV. See Fig. 5. 

Due to this heavy reliance on perturbative QCD, and consequently the rel- 
atively light reliance on experimental data, the uncertainty in Aa[ 5 j (m|) 
is reduced. The value quoted in Ref. (9) is 
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*<*hL m .(™l) = 0.02739 ± 0.00042, (17) 
which corresponds to 

a _1 (m z ) = 128.98 ± 0.06. (18) 
This disagrees with both Eq. 12 and 14 by about 0.10, but agrees with Eq. 16. 



The analysis of Eidelman and Jegerlehner 



The last analysis I will discuss is by Eidelman and Jegerlehner in Ref. (10). 
This work is another update of Ref. (7), and again trapezoidal integration is 
used. 

However, in order to take experimental errors into account properly, the 
data on R(s) was first processed in the following way: First, for all values 
of s with data points, a value of R(s) and its error was assigned to all the 
experiments in the region by linear interpolation between the closest data 
points belonging to that experiment. Then, the average of all the assigned 
values of R(s) was taken weighted with the experimental errors. The resulting 
weighted average was used in the integration. 

This analysis also used a more comprehensive set of data than any of the 
other analyses. 

After subtracting out the top quark contribution from the value reported 
in Ref. (10), we obtain 

*<*hLons(™l) = 0.02804 ± 0.00065, (19) 

and 

a _1 (m z ) = 128.89 ± 0.09, (20) 
which agrees very well with the previous estimate (7). 



Discussion 



While the large disagreement in a~ 1 (mz) of about 0.20 between Refs. (7) 
and (8) was disturbing, the new analyses by both authors (10,11) have de- 
creased the disagreement to less than 0.10. With the estimate of Ref. (9) 
falling in the same ballpark, around 0.10 would be a good estimate of the 
actual uncertainty in a~ 1 (mz). The smaller errors quoted in Refs. (9) and 
(11) are most probably underestimates due to the stronger assumptions they 
make on the functional form of R(s). 

A conservative conclusion would be that the current experimental data on 
R(s) does not allow for a substantially better determination of a(mz) than 
that already given in Ref. (7). Any further improvement on the uncertainty in 
a(mz) requires better measurements of R(s) in the low energy region below 
10 GeV. 
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TABLE 1. Comparison of the reliance on perturbative QCD in the evaluation of 
a(mz) between different authors. 



Author 


Ref. 


Energy Range (GeV) 


Order in a s 


a s (m z ) 


Burkhardt et al. 


(5,6) 


40 


— oo 


2 


0.12 ±0.02 


Jegerlehner 


(7) 


40 


— oo 


2 


0.117 ±0.010 


Swartz 


(8,11) 


15 


— oo 


3 


0.125 ±0.005 


Martin fc Zeppenfeld 


(9) 


3 - 3.9 


6.5 — oo 


3 


0.118 ±0.007 


Eidelman fc Jegerlehner 


(10) 


40 


— oo 


3 


0.126 ±0.005 



R 3 




MARK I (x 0.83) } 

Crystal Ball '90 \ 
(x 1.06) 



i ■ ■ 1 ■ i 

V(3.77) 



(b) 



R 3 



Crystal Ball '86 } 
Crystal Ball '90 <f 
(x 1.06) 



4 5 6 7 

Vs [GeV] 




I ' ' ' ' I ' ' ' ' 
*(3.77) (d) 



PLUTO (x 0.96) { 



3 4 5 6 

Vs [GeV] 



FIG. 5. The rescaling of data in Ref. (9). 
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TABLE 2. Current status of the determination of a s (mz). 



Method of Measurement 


Reference 




TITC 


(23) 


U.ll 6 ± U.UUo 


T decay 


(24) 


n i no _i_n nm 
U.lUo ± u.ulu 


T i-i. * 

Lattice 






charmonium spectrum 


(26) 


U.lUo it U.UUD 


I spectrum 


(27) 


0.H5 ± 0.002 


LEP 






Z lineshape data 


(1) 


0.125 ±0.004 


Ri only 


(1) 


0.128 ±0.005 


event shapes and jet rates 


(2) 


0.123 ±0.006 


SLD 






event shapes and jet rates 


(28) 


0.120± 0.008 



THE VALUE OF a s (m z ) 

Due to the limited number of pages I have been allocated, I will go through 
this section briefly. 

Current status 

In the case of a s (mz), there exist more than one way to determine it's 
value. Ref. (16) gives a comprehensive overview of the many ways to measure 
a s (mz). In table 2, I list some of the most recent determinations of a s (mz) 
using various techniques. 

Since the experimental error on is smaller than the theoretical error on 
the prediction based on Eq. 1, we can assume the SM and use the value of 
to determine a s (mz) to a better accuracy than Eq. 1. This gives the value 
oi s {mz) = 0.128 ± 0.005 listed in table 2. A global fit to all the Z line shape 
data, including R^, gives a s (mz) = 0.125 ± 0.004, also listed in table 2. 

In this approach where one assumes the SM and fits the value of a s (mz) 
to the data, any new physics beyond the SM will manifest itself not as a dis- 
agreement between the experimental measurement and theoretical prediction 
of some observable, but as disagreements between different determinations of 
a s (m z ). 

Such a disagreement may already have been seen. As you can see from ta- 
ble 2, the 'low-energy' measurements such as Deep Inelastic Scattering (DIS) 
and T decay favor a value of a s (mz) close to 0.11 while the 'high energy' 
measurements at LEP favor a value around 0.125. While this discrepancy 
may not seem too significant, considering the fact that unlike the a(mz) case 
these are all independent determinations based on different experiments, it 
is nevertheless puzzling that the LEP values of a s (mz) are systematically 
higher than the low energy values. 

Shifman argues in Ref. (30) that a value of a s (mz) as high as 0.125 would 
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correspond to a QCD scale of about 500MeV which is disfavored from the 
success of QCD sum rules and that the value 0.11 corresponding to a QCD 
scale of about 200MeV is more likely. 

An additional support for the lower value of a s (mz) is the recent measure- 
ment of the ratio Ri, = T(Z — > bb)/T(Z — > hadrons) at LEP which was found 
to be 2.4(7 larger than the prediction for the SM (1). If this is a signal that 
T(Z — > bb) is larger than its SM value, then one should also see a disagreement 
between theory and experiment in = T(Z — > hadrons) /T(Z — > But 
since one does not see a noticeable difference when comparing Eqs. 2 and 3, 
a larger value of T(Z — > bb) is not favored by R^ and one is led to conclude 
that the deviation in Ri, is just a statistical fluctuation. 

However, if the value of a s (mz) were as small as 0.11, then the story would 
be different. For a s (mz) = 0.110, the SM prediction of R^ will be lowered 
to 20.679 and the difference from Eq. 3 will be just right to accommodate an 
excess in T(Z — > bb) implied by Ri,. In this case, the disagreement between 
the SM and experiment will be at the 4<7 level (29-32). 

SUMMARY 

Though several authors have attempted to determination a(mz) to a better 
accuracy than Eq. 1, any substantial improvement requires better R(s) data 
at low energies. 

For the a s (mz) case, the disagreement between the low and high energy 
determination could be a signal for new physics which also causes Ri, to deviate 
from the SM. 
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